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Abstract 



Van der Waals coefficients for the heteronuclear alkali-metal dimers of Li, 
Na, K, Rb, Cs, and Fr are calculated using relativistic ah initio methods aug- 
mented by high-precision experimental data. We argue that the uncertainties 
in the coefficients are unlikely to exceed about 1%. 
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Considerable attention has been given to the determination of the coefficients of the 
leading term of the van der Waals attractions of two alkali metal atoms because of their 
importance in the simulation, prediction, and interpretation of experiments on cold atom 
collisions, photoassociation and fluorescence spectroscopy There is strong interest in 

heteronuclear molecules formed by pairs of different alkali metal atoms. Experiments have 
been carried out on trap loss in mixtures of Na with K Rb 0,0, and Cs |]lT| and on 



optical collisions [0 in a Na-Cs mixture and on molecular formation |TB[. The mixtures 
of magnetically-trapped alkali metal atoms Na-Cs and Na-K have been proposed as a 



means to search for evidence of an electric dipole moment to test for violation of parity 



and time reversal symmetry. We extend here previous studies of the van der Waals 
coefficient between pairs of identical ground state alkali metal atoms to unlike ground state 
atoms. 

The leading team of the van der Waals interaction is given at an atom separation R 

by 

v^'iR) = (1) 

where Cq^ is the van der Waals coefficient. We use atomic units throughout. 
The van der Waals coefficient may be expressed as 

2 \{vA\DA\sA)?\{vB\DB\tB)\' 

' 3^ {E^ - E^^) + (E,^ - EB) ' ^ ^ 

where \va) is the ground state atomic wave function of atom A with energy E^, and similarly 
for atom B, and \sa) and Its) represent complete sets of intermediate atomic states with, 
respectively, energies E^ and Ej^ . The electric dipole operators are Da = Z^i^i^^, where 
rf is the position vector of electron i measured from nucleus A, Na is the total number of 
atomic electrons for atom A, and similarly for atom B. 

At this point the two-center molecular-structure problem is reduced to the determination 
of atomic matrix elements and energies. The dependence on one-center atomic properties 
becomes explicit when Eq. (^) is cast into the Casimir-Polder form 

3 r°° 

Cq^ = — / aA{ioj)aB{ioj)duj , (3) 

TT JO 

where aA{i^) is the dynamic polarizability of imaginary argument for atom A given by 

, 2^[e^-E^)\{v^\Da\s^)? 
aAiiijj) = - > : — - — , (4) 

and a(to' = 0) is the ground-state static dipole polarizability. In the limit of infinite frequency 
the function aAi}'^^ satisfies 

(■ \ «\ 

(Xa\%^) ' (5) 

as a consequence of the nonrelativistic Thomas-Reiche-Kuhn sum rule. 
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Modern all-order many-body methods are capable of predicting electric-dipole matrix 
elements for principal transitions and energies in alkali-metals to within errors approaching 
0.1% |]T8[. Many-body methods augmented by high-precision experimental data for principal 



transitions, similar to those employed in PNC calculations |jT9|, have led to a high-precision 
evaluation of dynamic dipole polarizabilities for alkali-metal atoms The values of Ce 



previously calculated for homonuclear dimers are in excellent agreement with analyses 
of cold-atom scattering of Na |2^, Rb 0, and Cs [§,0]. Here we employ the same methods 



to compute the van der Waals coefficients for heteronuclear alkali-metal dimers. 

Precise nonrelativistic variational calculations of Cg for Li2 have been carried out 
They provide a critical test of our procedures. We separate the dynamic polarizability 
into valence and core contributions, which correspond respectively to valence-electron and 
core-electron excited intermediate states in the sum, Eq. (^). In our calculations for Li we 
employ high-precision experimental values for the principal transition 2s — 2pj, all-order 
many-body data and experimental energies for 3pj and 4pj intermediate states, and Dirac- 
Hartree-Fock values for higher valence-electron excitations. The high-precision all-order 
calculations were performed using the relativistic linearized coupled-cluster method trun- 
cated at single and double excitations from a reference determinant p8| , |23[| . Contributions 
of valence-excited states above 4pj were obtained by a direct summation over a relativistic 
B-spline basis set [2|] obtained in the "frozen-core" {V^~^) Dirac-Hartree-Fock potential. 
Core excitations were treated with a highly-accurate relativistic configuration-interaction 
method applied to the two-electron Li"^ ion. For the heavier alkali-metals [|15[ the random- 
phase approximation pSf was used to calculate this contribution. 

The principal transition 2s — 2pj accounts for 99% of the static polarizability and 96% 
of the Li2 dispersion coefficient. In accurate experiments McAlexander et al. [^| reported 
a lifetime of the 2p state of 27.102(9) ns (an accuracy of 0.03%) and Martin et al. 
reported 27.13(2) ns. In our calculations we employ the more precise value from Ref. ^6); in 
the subsequent error analysis we arbitrarily assigned an error bar of twice the quoted value 
of Ref. p6[, so that the two experiments are consistent. 



The dynamic core polarizability of Li was obtained in the framework of the relativistic 
configuration-interaction (CI) method for helium-like systems. This CI setup is described 
by Johnson and Cheng ||2^, who used it to calculate precise relativistic static dipole po- 
larizabilities. We extended their method to calculate the dynamic polarizability a{iuj) for 
two-electron systems. The numerical accuracy was monitored by comparison with results of 
Ref. 1^ for the static polarizability of Li"*" and with the sum rule, Eq. (^, in the limit of 
large frequencies. Core-excited states contribute only 0.5% to Cg and 0.1% to a(0) for Li. 
Their contribution becomes much larger for heavier alkali metals. 

We calculated static and dynamic polarizabilities and used quadrature, Eq. (|^), to obtain 
the dispersion coefficient. The results are Cq = 1390 and a(0) = 164.0. There are two major 
sources of uncertainties in the final value of Cq — experimental error in the dipole matrix 
elements of the principal transition, and theoretical error related to higher valence-electron 
excitations. The former results in a uncertainty of 0.12%, and the latter much less. The 
result Cq = 1390(2) is in good agreement with the nonrelativistic variational result of Yan 
et al. Cq = 1393.39. The slight discrepancy between the two values may arise because 
in our formulation, the correlations of core-excited states with the valence electron were 
disregarded as were intermediate states containing simultaneous excitation of the valence 
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electron with one or both core electrons. On the other hand, Ref. did not account for 
relativistic corrections. Relativistic contractions lead to a smaller value of Cg and to better 
agreement between the present result and that of Ref. [^]. Similar error analysis for the 
static polarizability of Li leads to a(0) = 164.0(1), which agrees with the numerically precise 
nonrelativistic result of 164.111 ||2^. An extensive comparison with other published data for 
the values of a(0) and Cq for lithium is given in Ref. ||2^. For the heavier alkali metal atoms 
we followed the procedures of Ref. ITH] to calculate a{iu!) 



The results for Cs are illustrated 
in Fig. |l|. They indicate that while most of the contribution to Cq comes from the resonant 
transition at ~ 0.05 a.u. the core excitations are significant. 

Results and Conclusions — We evaluated the dispersion coefficients for various heteronu- 
clear alkali- metal dimers with the quadrature Eq. @. The calculated values are presented in 
Table |. Most of the contributions to C^^ come from the principal transitions of each atom. 
An analysis of the dispersion coefficient of unlike atoms yields the approximate formula 



r^AAr^BB 



AEa + AEb 
VAEaAEb 



(6) 



where the energy separations of the principal transitions are designated as AEa and AEb- 
Eq. (P) combined with the high-accuracy values of Cq for homonuclear dimers gives 
accurate approximations to our results based on Eq. (j^). For example, Eq. overestimates 
our accurate value from Table | for Li-Na by 0.4% and for Cs-Li by 2%. We may use Eq. (H) 
to estimate the uncertainties SC^^ in the heteronuclear cases from the uncertainties 5Cq^ 
and SCq^ in the homonuclear dispersion coefficients. 



'5C, 



AA^ 



+ 



1/2 



The accuracy of Cg for homonuclear dimers was assessed in Ref. 
sion for the Rb dimer is given in Ref, 



15 and a detailed discus- 
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Analyzing the error in this manner using the 
quoted coefficients and their uncertainties from Ref. |]15| we find that most of the disper- 
sion coefficients reported here have an estimated uncertainty below 1%. The corresponding 
values are given in parentheses in Table |I[ 

In Fig. 1^ we present for the dispersion coefficients of the dimers involving Cs a com- 
parison between our calculated values and the most recent determinations [§,0. We give 
the percentage deviation from our calculations. It is apparent that the other calculations 
that employed one-electron model potentials and accordingly omitted contributions from 
core-excited states yield values systematically smaller than ours. 

The discrepancies are most significant for Cs2 where the number of electrons is greatest. 
Fig. ^ also compares the values for the Cs2 dimer with values deduced from ultracold-collision 
data [|6|,pO[. The agreement of our prediction 6851(74) [|15| with their values for in Cs2 is 
close. Core-excited states contribute 15% l3l||l^ to the value of the Cg coefficient for the Cs 
dimer and are needed to fulfill the oscillator strength sum rule, Eq. (^. In the present ap- 
proach the contributions of core-excited states to dynamic polarizabilities are obtained using 
the random-phase approximation, which nonrelativistically satisfies the oscillator strength 
25[] . In the inset of Fig. 0, it is illustrated that our calculated aiioj) 
as LO becomes asymptotically large, where = 55 for Cs. While the 



sum rule exactly 
approaches NjuP' 
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deviation between the present calculations and the model potential calculations are smaller 
for dimers involving lighter atoms, an accurate accounting of core-excited states is essential 
to achieve high accuracy in dispersion coefficient calculations for heavy atoms [pl]-p3[]. 

Few experimental data are available for comparison in the heteronuclear case, except for 
NaK. The results from investigations of NaK molecular potentials based on spectral analy- 
sis 0] are compared to our value in Table |I[ Our value is smaller than the experimental 
values. Earlier theoretical calculations of dispersion coefficients for NaK have been tabu- 
lated and evaluated by Marinescu and Sadeghpour ||35| and by Zemke and Stwalley |^ 



Those values are generally lower than our value of 2447(6) except for that of Maeder and 
Kultzelnigg [0 who give 2443. 

The present study extends the application of modern relativistic atomic structure meth- 
ods to calculations of ground state van der Waals coefficients of Li2 and of the heteronuclear 
alkali-metal atoms. We argue that the uncertainty of the coefficients is unlikely to exceed 
1%. Additional experimental data from future cold-collision experiments or spectroscopy 
would provide further tests of the present calculations. 

This work was supported by the Chemical Sciences, Geosciences and Biosciences Division 
of the Office of Basic Energy Sciences, Office of Science, U.S. Department of Energy and by 
the National Science Foundation under grant PHY97-24713. The Institute for Theoretical 
Atomic and Molecular Physics is supported by a grant from the NSF to Harvard University 
and the Smithsonian Institution. 
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FIG. 1. The dependence of the dynamic dipole polarizabihty a{iio) with frequency to for Cs. 
The inset illustrates the behavior of the quantity uj'^a{iLo) at asymptotically large lo, where the 
dashed line represents the contribution of the core-excited states to the total w^a(iu;) (solid line) 
and the arrow marks the non-relativistic limit = 55 following from the sum rule, Eq. (|). All 
quantities are in atomic units. 
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LiCs NaCs KCs RbCs Csj 
FIG. 2. Percentage deviation of results of recent calculations |35,37| from our values for van 
der Waals coefficients Cq for Cs-Li, Cs-Na, Cs-K,Cs-Rb, and Cs-Cs. The values with error bars 
placed along the horizontal line at correspond to our results with the estimated uncertainties. 
Circles represent the results of Ref. and triangles the results of Ref. [37|. For Cs-Cs, to the 
right of the vertical dotted line, we show the difference between our present prediction, our earlier 
prediction |15| and the values deduced from cold-collision data in Ref. (square) and Ref. |3C] 
(diamond) . 
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TABLES 



TABLE L Dispersion coefficients Cq and their estimated uncertainties (parentheses) for aL 



kaU-metal 
Ref. 


atom pairs in atomic units. 


Coefficients 


for Na2, K2, 


Rb2, Cs2, and Fr2 


are from 




Li Na 


K 


Rb 


Cs 


Fr 


Li 


1389(2) 1467(2) 


2322(5) 


2545(7) 


3065(16) 


2682(23) 


Na 


1556(4) 


2447(6) 


2683(7) 


3227(18) 


2842(24) 


K 




3897(15) 


4274(13) 


5159(30) 


4500(39) 


Rb 






4691(23) 


5663(34) 


4946(44) 


Cs 








6851(74) 


5968(60) 


Fr 










5256(89) 
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TABLE II. Comparision of present theoretical and experimental values for the dispersion 
coefficient for NaK. 



Reference 



This work 

Russier-Antoine et al, [34| 
Ishikawa et al, [35] 
Ross et al, [39| 



2447(6) 
2519(10)^ 
2646(31)'^ 
2669.4(20)'^ 



^Experiment. 



9 



REFERENCES 



H. M. J. M. Boesten, C. C. Tsai, J. R. Gardner, D. J. Heinzen, and B. J. Verhaar, Phys. 
Rev. A 55, 636 (1997). 

J. L. Roberts, N. R. Claussen, J. P. Burke Jr., C. H. Greene, E. A. Cornell, and C. E. 
Wieman, Phys. Rev. Lett. 81, 5109 (1998). 

C. J. Williams, E. Tiesinga, P. S. Julienne, H. Wang, W. C. Stwalley, and P. L. Gould, 

Phys. Rev. A 60, 4427 (1999). 

J. L. Bohn, Phys. Rev. A. 61, 053409 (2000). 

C. Amiot and J. Verges, J. Chem. Phys. 112, 7068 (2000). 

P. J. Leo, C. J. Williams, and P. S. Julienne, Phys. Rev. Lett. 85, 2721 (2000). 

M. S. Santos, P. Nussenzveig, L. G. Marcassa, K. Helmerson, J. Flemming, S. C. Zilio, 

and V. S. Bagnato, Phys. Rev. A 52, R4340 (1995). 

M. S. Santos, P. Nussenzveig, A. Antunes, P. S. P. Cardona, and V. S. Bagnato, Phys. 
Rev. A 60, 3892 (1999). 

G. D. Telles, I. G. Marcassa, S. R. Muniz, S. G. Minanda, A. Antunes, C. Westbrook, 
and V. S. Bagnato, Phys. Rev. A 59, R23 (1999). 

Y. E. Young, R. Ejnisman, J. P. Shaffer, and N. P. Bigelow, Phys. Rev. A 62, 055403 
(2000). 

J. P. Shaffer, W. Chalupczak, and N. P. Bigelow, Phys. Rev. A 60, R3365 (1999). 
J. P. Shaffer, W. Chalupczak, and N. P. Bigelow, Phys. Rev. A 61, 011404 (1999). 
J. P. Shaffer, W. Chalupczak, and N. P. Bigelow, Phys. Rev. Lett. 82, 1124 (1999). 
R. Ejnisman, P. Rudy, N. P. Bigelow, P. S. P. Cardona, A. M. Tuboy, D. M. B. P. 
Milori, V. S. Bagnato, and I. D. Goldman, Braz. J. Phys. 27, 247 (1997). 
A. Derevianko, W. R. Johnson, M. S. Safronova, and J. F. Babb, Phys. Rev. Lett. 82, 
3589 (1999). 

J. Goodisman, Diatomic Interaction Potential Theory (Academic Press, New York, 
1973), Vol. 2. 

A. Dalgarno, in Intermolecular Forces, Vol. 12 of Advances in Chemical Physics, edited 
by J. O. Hirschfelder (Wiley, New York, 1967), p. 143. 

M. S. Safronova, W. R. Johnson, and A. Derevianko, Phys. Rev. A 60, 4476 (1999). 

S. A. Blundell, W. R. Johnson, and J. Sapirstein, Phys. Rev. Lett. 65, 1411 (1990). 

F. A. van Abeelen and B. J. Verhaar, Phys. Rev. A 59, 578 (1999). 

C. Chin, V. Vuletic, A. J. Kerman, and S. Chu, Phys. Rev. Lett. 85, 2117 (2000). 

Z.-C. Yan, J. F. Babb, A. Dalgarno, and G. W. F. Drake, Phys. Rev. A 54, 2824 (1996). 

S. A. Blundell, W. R. Johnson, Z. W. Liu, and J. Sapirstein, Phys. Rev. A 40, 2233 

(1989). 

W. Johnson, S. Blundell, and J. Sapirstein, Phys. Rev. A 37, 307 (1988). 

M. Amusia and N. Cherepkov, Case Studies in Atomic Physics 5, 47 (1975). 

W. I. McAlexander, E. R. I. Abraham, and R. G. Hulet, Phys. Rev. A 54, R5 (1996). 

F. Martin, M. Aubert-Frecon, R. Bacis, P. Crozet, C. Linton, S. Magnier, A. J. Ross, 

and I. Russier, Phys. Rev. A 55, 3458 (1997). 

W. R. Johnson and K. T. Cheng, Phys. Rev. A 53, 1375 (1996). 

M. S. Safronova and W. R. Johnson, available at 

|http:/ /www.nd.edu/~johnson/Publications/C6_Rb.pdl| (unpublished). 



10 



[30] C. Drag, B. L. Tolra, B. T'Jampcns, D. Comparat, M. AUegrini, A. Crubellier, and P. 

Fillet, Phys. Rev. Lett. 85, 1408 (2000). 
[31] A. Dalgarno and W. D. Davison, Mol. Phys. 13, 479 (1967). 
[32] F. Maeder and W. Kutzelnigg, Chem. Phys. 42, 95 (1979). 

[33] M. Marinescu, H. R. Sadeghpour, and A. Dalgarno, Phys. Rev. A 49, 982 (1994). 
[34] I. Russier-Antoine, A. Ross, M. Aubert-Frecon, F. Martin, and P. Crozet, J. Phys. B 
33, 2753 (2000). 

[35] M. Marinescu and H. Sadeghpour, Phys. Rev. A 59, 390 (1999). 

[36] W. T. Zemke and W. C. Stwalley, J. Chem. Phys. Ill, 4956 (1999). 

[37] S. H. Patil and K. T. Tang, J. Chem. Phys. 106, 2298 (1997). 

[38] K. Ishikawa, N. Mukai, and W. Tanimura, J. Chem. Phys. 101, 876 (1994). 

[39] A. J. Ross, C. Effantin, J. d'Incan, and R. F. Barrow, Mol. Phys. 56, 903 (1985). 



11 



